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ABSTRACT 

We have calculated the electric dipole strength distributions in the unstable neutron 
rich oxygen isotopes 18 > 20 > 22 O, in a model which include up to four quasi-particle-type 
configurations. The model is the extension, to include the effect of the pairing corre- 
lations, of a previous model very successful around closed shell nuclei, and it is based 
on the quasi-particle-phonon coupling. Low-lying dipole strength is found, which ex- 
hausts between 5 and 10% of the Thomas-Reiche-Kuhn (TRK) energy-weighted-sum- 
rule (EWSR) below 15 MeV excitation energy, in rather good agreement with recent 
experimental data. The role of the phonon coupling is shown to be crucial in order to 
obtain this result. 
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1 Introduction 



The study of Giant Resonances (GR) in nuclei far from the stability line, such as halo 
and skin nuclei 0, is still in its infancy. The observation of large El strength at 
low excitation energy, just above threshold in light drip-line nuclei, has produced strong 
excitement and large amount of activities, both experimental and theoretical. Electro- 
magnetic dissociation measurements of this low-lying stength have been performed for 
the neutron halo nuclei 6 He [§, 8 He U, u Li §, n Be [|, 12 Be§, 14 Be[|, 19 C @, and for 
the proton halo nuclei 8 B |T(J and 13 From the theory side, a quantum-mechanical 
threshold effect was shown to enhance the transition strength from the loosely bound 



nucleons to the continuum |T2|. [Tj|], essentially produced by the optimal matching of the 
wavelength of the scattering states with the large penetration length into the forbidden 
region of the weakly bound orbitals. In the case of two-nucleon halo nuclei, the en- 
hancement due to the coherence in the transition amplitudes between the loosely bound 
nucleons, and also between these and the core, has been found important in recent shell- 
model calculations with the use of extended wave functions for the matrix elements of 
the dipole operator [jrj], |TJ] (see also [[To]). 

Moving to slightly heavier, eventually skin nuclei, the first attempt [|17|] of a system- 
atic experimental study of GR has been recently performed at GSI, with the investi- 
gation of the dipole response of unstable neutron rich oxygen isotopes up to 22 and 
to an excitation energy of 30 MeV. Use has been done of the electromagnetic excita- 
tion in heavy-ion collisions at around 600 MeV/u beam energy. Low-energy collective 
modes in 20 O were also measured recently for the first time by means of inelastic proton 



scattering [|TJ . 

Mean time, theoretical studies have been performed in the context of mean-field theo- 
ries, like self-consistent Hartree-Fock and Random Phase Approximation (HF+RPA) cal- 
culations fl9| , eventually including exactly the crucial coupling to the continuum p0[ , 
or HF-BCS plus quasi-particle RPA (QRPA) calculations p2[ (which include the contri- 
bution of the pairing force), as well as relativistic RPA (RRPA) calculations |23||, in order 



to investigate how the strength distributions of different multipole operators, in partic- 
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ular monopole, dipole and quadrupole, are evolving as one approaches the drip lines. 
The main features result to be: (a) the unusual concentration of multipole strength at 
the continuum threshold, already mentioned, when the threshold energy becomes of the 
order of one or few MeV, and (b) the strong interplay or mixing between the so-called 
isoscalar and isovector modes familiar from the response of stable nuclei. 

Much less, and not at all systematically, has been done beyond mean field, that is, 
including the coupling of single-particle states to vibrational modes, in particular low- 
lying vibrations |[24|| , as widely done in the case of stable nuclei in the last decades. Some 
pioneering works on the effects of this coupling on the single-particle self-energy and the 
GR strength function in netron-rich nuclei are reported, e.g., in Refs. [^5, 26, 27, 28, 2E], 
Of course, these correlations beyond mean field are also included, although somewhat 
hidden, in the shell- model calculations performed by several groups In particular, 
we mention the calculation of Ref. where the dipole strength functions in the O 
isotopes were calculated, before the results of the GSI experiment quoted above were 
known. 

Thus, it appears timely to extend in the present paper the microscopic model de- 



scribed in detail in Ref. |32| and aimed to a microscopic description of the excitation and 
decay of GR in stable nuclei, to the case of the neutron rich systems. This is done in the 
present work, were pairing correlations are included in the model in a simple way. The 
model is based on HF-RPA with Skyrme effective forces and includes on top escape and 
spreading effects in a consistent way. It was found to give accurate predictions, around 
closed shell nuclei, for the excitation energy and strength of GR as well as for widths 



and particle decay branching ratios p3 |. 

The paper is organized as follows. In Section 2, the model and its extension are 
presented. The application to the calculation of the dipole strength functions in the O 
isotopes is described in Sect. 3, including also the comparison with the recent experi- 
mental data of Ref. |T7| and with the shell-model results of Ref. J31|. Preliminary results 
were published in conference proceedings ||34j| . Summary and conclusions are drawn in 
Sect. 4. 
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2 Description of the model 



In this section, we shortly present the extension of the model of Ref . , to include in a 
minimal way the effects of the pairing interaction. Indeed, we are interested to perform 
calculations along isotopes chains and the inclusion of pairing correlations is called for. 
The starting input is an effective Hamiltonian H which includes a two-body interac- 



tion of the Skyrme type |35[ . The self-consistent mean field is determined by solving the 
Hartree-Fock (HF) equations on a radial mesh, coupled with the standard BCS equations 
for the added pairing interaction (within the constant pairing gap approximation). This 
method is well known for spherical nuclei in the case of Skyrme interactions [[35], j36fl . 



The unoccupied states, at negative as well as positive energy, are determined by using 
a harmonic oscillator basis, so the continuum is discretized. This procedure defines a 
subspace labelled by Qi and including discrete two quasi-particle (2qp) configurations. 
The standard |37| QRPA equation in matrix form are solved on a basis of 7^=1" con- 
figurations belonging to Qi, to obtain collective (or non-collective) low- lying and GR 
modes (see also Ref. fl2"2|| ). Special attention is paid to project out of the spectrum the 



spurious center-or-mass state, by renormalizing by a few percent the matrix elements of 
the residual interaction in such a way that the spurious state is pushed to zero energy. 

To account for the spreading width we build a second subspace Q 2 (the subspace 
P introduced in Ref. [E| to account for the escape width is not included in the present 



model, see also next section). In Q 2: there are the main configurations known ]23] to play 
a major role in the damping process of GR: they are 2qp states coupled to a collective 
vibration calculated using the discrete QRPA in the Qi subspace as described above. 
We discuss in the next section what vibrations have been selected and included in Q 2 . 

Using the projection operator formalism, it is easily found that the effects of coupling 
the subspace Q 2 to Q\ are described by the following effective Hamiltonian H, acting in 
the Qi space 



H{E) = Q l HQ 1 + W i {E) 

QiHQi + QiffQ 2 1 . Q 2 HQ U 

E - Q 2 HQ 2 + ie 

(1) 
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where E is the excitation energy. For each value of E, the QRPA equations corresponding 
to this complex Hamiltonian 7~C{E) are solved (this is done using the basis of the discrete 
QRPA states, eliminating those with negligible strength). The resulting sets of complex 
eigenstates \v) and eigenvalues E v — iT u /2, enable to calculate all relevant quantities, in 
particular the strength function S(E) of the operator F of interest: 

More in detail, each matrix element of W^(E) written in the basis of the discrete 
QRPA states is a linear combination of matrix elements W^yg(E) in the basis of the 
2qp configurations (a/3), (jd) etc. These latter matrix elements are sums of the eight 
terms depicted in Fig. 1. To evaluate these diagrams, we employ the following expression 
for the particle-vibration coupling Hamiltonian, 

V = Y. E (a\el L) (r)v(r)Y LM (r)\(3) a^ap, (3) 

a/3 LnM 

and we use the standard BCS expressions to relate the quasi-particle operators c a and c a 
to the single-particle ones a£, and a a (see Eq. fljjD ). In Eq. (||) the vibration (phonon) |n) 
is characterized by its angular momentum L and its transition density Q^(r), while the 
form factor v(r) is related to the particle- hole interaction V p h derived from the Skyrme 
force as V p h{fi,r 2 ) = v(rx)S( , rx, r 2 ). The detailed expressions of the diagrams are given 
in the Appendix. 

3 Application to the dipole response of 1822 

We use the model discussed above to calculate the electric dipole strength distributions 
in the unstable neutron rich oxygen isotopes. As anticipated in the introduction, in a 



recent experiment [17j at GSI its knowledge has been extended to the isotopes heavier 



than 18 O, for which it was well known |38|, with systematic measurements from 



^4=17 to A=22. Low-energy strength was detected, which exhausts up to 12% of the 
classical Thomas-Reiche-Kuhn (TRK) energy-weighted sum rule (EWSR) at excitation 
energies below 15 MeV (approximate threshold for the emission of the protons, which 
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are not detected in the experiment). The comparison with these data constitutes an 
adequate test for our model. 

The calculations have been performed in the isotopes 18 > 20 > 22 O, and using the effective 
force SIII [ p5|| . This force, in the HF approximation, reproduces very well the binding 
energy of these nuclei (although it fails to fix the drip line for Z=8 at A=2A). However, 
the corresponding self-consistent RPA description of the low-lying vibrations is rather 
poor in the comparison with the experimental data, and is improved by the inclusion of 



some pairing correlations |1| ^2] . Therefore, we do include a pairing force with constant 
matrix elements of strength G = 0.4 MeV (about 30% of the standard value), and adopt 
the corresponding values of A n for the neutrons in our constant gap calculation. These 
pairing gaps are smaller than 1 MeV (0.74, 0.79 and 0.21 MeV respectively in 18 > 20 > 22 O), 
and the binding energies differ only by about 1% from the HF value. The B(E2) values 
for the low-lying 2 + states are somewhat improved compared to the RPA calculation 
(in the case of 20 O, the experimental value for the B(E2) is 28±2 e 2 fm 4 [10]] , RPA gives 
only about 3 e 2 fm 4 and our QRPA with the small gap gives 15.6 e 2 fm 4 ). 

The HF-BCS calculations are done in coordinate space, using a radial grid of 0.1 fm 
extending up to 15 fm. The use of a constant pairing gap A leads to unrealistic results 
unless a cutoff is introduced in the single-particle space, in such a way that states above 
this cutoff do not feel any pairing interaction: in our case, the cutoff is just above the 
ld3/2 neutron state. In the QRPA calculations for the dipole and for the isoscalar modes 
to be included in the subspace Q2, the residual interaction between quasi-particles is 
derived from the Skyrme force without including the pairing contribution. The 2qp 
basis is chosen large enough so that the appropriate EWSR are satisfied. In the case of 
the dipole, we refer to the double commutator sum rule which is enhanced with respect 
to the classical TRK sum rule by a factor 1.33 in the case of the SIII interaction at hand. 

In the QRPA calculation, the dipole strength below 15 MeV is found smaller than the 
experimental one (see Table p]). Then, the full coupling to the complex configurations of 
the Q2 space is taken into account. The phonons included are those of multipolarity 1~, 
2 + , 3~ and 4 + which lie below 30 MeV and absorb more than 5% of the total strength. 
(In the case of 18 O, the + modes have been included but their effect is found to be 
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rather small). The final results for the integrated cross section up to 15 MeV is shown 
in Table |l| in comparison with the experimental data and with the results of the 



shell-model calculation of f31fl . The coupling with phonons increases the low-energy 
cross section, bringing it in rather good agreement with the data. Moreover, the drop 
of the cross section in 22 O is also reproduced. In the case of 18 O our results are also 
remarkably similar to the shell-model ones, suggesting that we are including the most 
important correlations. 

The entire cross sections obtained with the full coupling for the three isotopes are 
displayed in Fig. 2. The progressive cumulation of cross section at lower energies with 
increasing A is evident, as well as the mild dependence on A of the main GDR peak and 
the appearance of more fragmented and peculiar structures at higher energies. In the case 
of 18 0, the peak energy at 24.2 MeV compares very well with the known experimental 
value |5Bfl , although the full experimental strength distribution looks wider (see Fig. 3(e) 

of H). 

These features are not qualitatively very different in the calculations we performed 
with much larger pairing gap values, that is, the standard 12 / y/A MeV values. The 
results of these latter calculations are displayed in Fig. 3. In this case, compared to 
the calculation with the small gaps, the GDR main peak is more damped and the high 
energy structures more pronounced. The largest difference appears for the case of 22 0: 
moreover, in this nucleus the strength below 15 MeV is definitely larger than in the 
small gap case, at variance with 18 ' 20 O where this strength is not very sensitive to the 
value of A ra we have employed (see Table [I]). The fact that the GDR peak energy is not 
markedly affected by the size of the pairing gap in all cases can be understood since this 
is essentially ruled by mean field effects and at that level, the effect of pairing should not 
be important for states at energies much larger than « 2 A. Concerning the coupling 
with the phonons, for large values of the pairing gap, the low-energy surface vibrations 
become more collective but tend also to be pushed at higher energies (see also Fig. 1 
of Ref. P2| ) and these two facts have opposite consequences on the effectiveness of the 
phonon coupling. In the present case we understand from Fig. 3 that the first of these 
two effects (increased phonon collectivity) seems to be somewhat dominant. 
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We have not discussed so far the sensitivity of the results to the choice of the effective 
Skyrme force. An example can be seen in Fig. where the full and dashed lines are two 



calculations for with the large gap, using respectively the forces SIII and SGII [41 



The difference between the two curves is essentially accounted for by a downward shift 
of about 2 MeV. This is mainly due to the different effective mass associated with the 
two interactions. 

The most important outcome of our study is that the effect of the coupling with 
complex configurations including one phonon is very large and unavoidable if the com- 
parison with experimental data is envisaged. The point is clearly emphasized in Fig. 5, 
for 18 0. The mean field (i.e., QRPA) result, with the simple smearing by means of 1 MeV 
width Lorentzian functions, are qualitatively altered when the coupling with phonons 
is included. This leads to a largely spread cross section, in much better agreement 
with the experiment. The effect is certainly larger than in stable, heavy nuclei, where 
most of our experience is. One main reason is the lack of the large cancellation pi | 
between the diagrams in the upper and lower two lines of Fig. 1. This is due to the very 
asymmetric phase-space available, in these light nuclei, for the hole-like and particle-like 
contributions to and was already discussed in Ref. [42] for the shell- model calculation 
of the giant quadrupole resonance (GQR) in 16 compared to 40 Ca. This asymmetry 
goes together, in hindering large cancellations, with the small number of components in 
the wave functions of the low-lying structures, as reported in Table and in Table [j| 
Consequently, these structures should not be considered of collective nature. 

Admittedly, this small number of components makes more critical than in the heavier 
nuclei the omission of the coupling to the continuum. The neutron separation energy in 
22 is still large, of the order of 7 MeV, but certainly the large role played in the wave 
functions by neutron quasi-particles like the 2p 3 / 2 in the continuum poses some questions, 
worthwhile to be investigated. In any case, we show in Fig. 6 the results at the level of 
HF (that is, the unperturbed dipole strength function) obtained either taking properly 
the continuum into account or discretizing it. By choosing reasonable parameters to 
average the discrete results (see the caption to the figure) we reproduce very well the 
exact results, and this means that peak regions and corresponding intensities can be, in 
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this case, accounted for by the discrete calculation. 



4 Conclusion 

In the present work we have calculated the electric dipole strength distributions in the 
unstable neutron rich oxygen isotopes 18 > 20 > 22 O. The model we used is the extension, to 



include the effect of the pairing correlations, of a previous model [32] which has been very 
successful around closed shell nuclei, and which is based on the quasi-particle-phonon 
coupling. Thus, the 2qp configurations of the standard QRPA approach, are coupled 
to 4qp-type configurations built with two uncorrelated quasi-particles plus a collective 
phonon (these are known to be particularly efficient in redistributing the GR strength). 

Low-lying dipole strength is found, exhausting between 5 and 10% of the TRK sum 
rule below 15 MeV excitation energy, in rather good agreement with the recent exper- 



imental data of ||1 7| . In order to obtain this result, and more generally to spread out 
the whole strength distribution including the main GDR peak, the role of the coupling 
with phonons appears to be crucial: therefore, the importance of particle-vibration cou- 
pling in neutron-rich nuclei is one of the most important outcomes of our paper. The 
comparison with the strength distributions obtained by means of large scale shell-model 
calculations f3]J is satisfactory in terms of integrated quantities, especially in 18 O, al- 
though the detail of the distributions looks quite different. 

Confident in our model, we are planning to apply it in the case of much heavier 
nuclei, relevant, e.g., for the study of the r-process and nucleosyntesis. In these nuclei, 
other types of microscopic calculations like the full shell-model may result prohibitive 
for computational reasons and, to describe the low-lying dipole strength, macroscopic 
models [[43], |44| which may not contain all proper physical ingredients have been used so 
far 
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Appendix 

In this Appendix we provide the expressions of the eight diagrams contributing to the 
matrix element s (E) mentioned in Sec. 2 and depicted in Fig. 1. 



W\ = 5(a, -y)5{jp,m + 5(a, /^-^[l + 5(a, 5)]- 1/2 (2 J/3 + l)" 1 £ [1 + 5(a, a)]" 1 / 2 



Wij = 8(a,8)5(j f3 ,j,)[l + 5(a,(3)]- 1 / 2 [l +6(a, 1 )]- 1 / 2 (2j p + l)" 1 £ [1 + 6(a,a)]-^ 



W\ n = WW«J 7 )[l + *(«.0)]^ 



W\ v = 5((3, -y)6(j a ,j,)[l + 5(a, (3)]- l / 2 [l + 5((3, d)]- 1 ^ + l)- 1 £ [1 + 5((3, a)]' 1 / 2 
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= [1 + 5(a, P)]- l/2 [l + 5( 7 , 5)]- 1 / 2 [l + 5(a, 5)}' 1 / 2 (U a U, - Kig^ - V P V S ) 
J2\ Ja 3P t f [ dru a (r)u^r)eip(r)v(r) [ drup(r)u 5 (r)g£\r)v(r) 

Ln { 3 5 h ^ J J J 

(-l)^ +J («||r L || 7 ) • (P\\Yl\\8) ■ - + ^ + c n ) + (8) 



(_l)^( a ||y L ||5) • </3||Xc|| T > • [u - (E a + Ep + u^+irj]- 1 ; (9) 



W^jj = [1 + <S(a,/3)]- 1/2 [l + 5( 7 ,5)]- 1/2 [l + 5(/5, 7 )]~ 1/2 (f/af/ 7 - VaVjiJJpUs - V P V 5 ) 
E| ^ J £ }/ dru Q (r)« 7 (r)^ L) (Ov(r) / drup(r)u s (r) Q ^ (r)v(r) 
(-l)w + *r +J (a||y L || 7 > • (/Vl||5> • [«- + E 7 + w„) (10) 



E{ j ^ I } Jdru^u^g^irMr) f dru a (r)u s (r)gl L \r)v(r) 
(-iy^(a\\Y L \\5) ■ ((3\\Yl\\i) • - (^ + ^ + ^n) (11) 

In the above expressions, a labels an intermediate quasi-particle state, E and u(r) 
indicate respectively the quasi-particle energy and the single particle radial wave function 
whereas U and V are the usual BCS coefficients which relate single particle and quasi- 
particle operators according to 

4 = U a al + V a a s , (12) 

where a is the time-reversed of a. Concerning the vibrational states (phonons), these are 
labelled by the multipolarity L and the index n: tu n and Sg^(r) are the corresponding 
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energies and radial transition densities as already discussed in the text in Sec. |2| (v (r) 
is also discussed there). The averaging parameter t] is set at 0.5 MeV in the present 
calculation. 

The above expressions are derived writing the single particle wave functions as 
{unii^) / r)i l Yi m {f) (phase convention II of Ref. H37] ). In the limit of zero pairing gap, 



the expressions reported in Ref. [[3^] are recovered, where however, phase convention I 
of Ref. f|37|l without the factor i l , was used. 
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Fig. 1. Diagrams corresponding to the coupling of the 2qp configurations to the more 
complicated states including a phonon. See Sec. 2 for the discussion of this coupling, 
the detailed analytic expressions being given in the Appendix. 

Fig. 2. Total photoabsorbtion cross section for the isotopic chain 18 > 20 > 22 O obtained 
within the full phonon coupling (using small pairing gaps as discussed in Sec. 3). 

Fig. 3. Same as Fig. 2 in the case of large pairing gaps (i.e., A n =12/VA MeV). 

Fig. 4. Total photoabsorbtion cross section in the case of 18 and of the large pairing 
gap, using two different effective forces namely SIII and SGII. 

Fig. 5. Total photoabsorbtion cross section for 18 0. The large spreading induced by the 
coupling with the phonons included in the complete calculation (full line) contrasts with 
the sharp peaks of the QRPA result (dashed line). The integrals of the two curves over 
the entire range are respectively 303 and 320 MeV-mb. 

Fig. 6. Unperturbed (i.e., Hartree-Fock) dipole strength obtained for 22 either using 
the proper continuum (full line) or averaging (dashed line) the sharp peaks (bars) of the 
discretized continuum. For the averaging procedure, Lorentzian functions having width 
2.2 MeV and 1.6 MeV have been employed, respectively for peaks below or above 12.5 
MeV. Units are fm 2 for the bars and MeV _1 -fm 2 for the lines. 



Table 1: Total photoabsorbtion cross section integrated up to 15 MeV obtained in our 
QRPA (second row and third row) and QRPA plus phonon coupling (third and fourth 
row) calculations. The different choices for the pairing gap A n (either small or large) are 
discussed extensively in the text. The results are compared with the photoneutron cross 
section measured at GSI (first row) and with the shell-model calculation of Ref. |JT| 
(last row). All numbers are in MeV-mb. 





18q 


2U 


22 


Experiment 








22.14 ± 2.4 


33.7 ± 3.5 


23.22 ± 2.75 


QRPA (small gap) 








11.43 


19.28 


15.07 


QRPA (large gap) 








11.00 


18.21 


22.36 


QRPA plus phonon 


coupling 


(small 


gap) 


18.26 


25.65 


17.05 


QRPA plus phonon 


coupling 


(large j 


?ap) 


18.55 


25.92 


22.87 


Shell model 








17.14 


31.26 


30.40 
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Table 2: Main components of the wave functions of the peaks around 12 MeV in 18 > 20 O 
(visible in Fig. 2). Only the real parts of the amplitudes which are larger in absolute value 
than 0.1, are reported, n and v label the proton and neutron amplitudes, respectively. 





2UQ 22Q 


Z/,lpi /2 2Si/ 2 
I>,ld 5 / 2 lf 7 /2 

z/,ld 5/2 2p 3/2 

7T,lpi/ 2 2Si/ 2 

z/,2si/ 2 2p 3/2 
z/,2si /2 2pi /2 


0.82 0.83 
0.31 0.35 

-0.32 

-0.18 

0.30 
0.10 



Table 3: Same as Table 1, in the case of the peak around 14.2 MeV in 20 O. 





2u Q 


Z/,ld 5/2 lf 7 /2 

z/,2si /2 2p 3/2 

Z/,lp 3/2 ld 5/2 


0.18 
-0.63 
-0.18 
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Large gaps 
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